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Abstract 

We examine the effects of thermal fluctuations on thin elastic filaments with non- 
circular cross-section and arbitrary spontaneous curvature and torsion. Analytical 
expressions for orientational correlation functions and for the persistence length of 
helices are derived, and it is found that this length varies non-monotonically with 
the strength of thermal fluctuations. In the weak fluctuation regime, the local helical 
structure is preserved and the statistical properties are dominated by long wavelength 
bending and torsion modes. As the amplitude of fluctuations is increased, the helix 
"melts" and all memory of intrinsic helical structure is lost. Spontaneous twist of 
the cross-section leads to resonant dependence of the persistence length on the twist 
rate. 

1. Introduction 

Modern polymer physics is based on the notion that while real polymers can be arbitrarily 
complicated objects, their universal features are captured by a minimal model in which poly- 
mers are described as continuous random walks. While this approach has been enormously 
successful and led to numerous triumphs such as the understanding of rubber elasticity!, the 
solution of the excluded volume problem and the theory of semi-dilute polymer solutions^, 
it is ill-suited for the description of non-universal features of polymers which may depend 
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on their chemical structure in a way that can not be captured by a simple redefinition of the 
effective monomer size or its second virial coefficient. For relatively simple synthetic poly- 
mers, such "local details" can be treated by polymer chemistry -type models (e.g., rotational 
isomer state modelB). However, chemically-detailed approaches become prohibitively diffi- 
cult (at least as far as analytical modeling is concerned) in the case of complex biomolecules 
such as DNA, proteins and their assemblies and a new type of minimal models is needed to 
model recent mechanical experiments on such systemsiH^. Such an alternative approach is 
to model polymers in the way one usually thinks of them, i.e., as continuous elastic strings 
or filaments which can be arbitrarily deformed and twisted. However, while the theory of 
elasticity of such objects is well-developed!!!, little is known about the statistical mechanics 
of fluctuating filaments with arbitrary natural shapes. The main difficulty is mathemati- 
cal in origin: the description of three-dimensional filaments with non-circular cross-section 
and non-vanishing spontaneous curvature and twist@, involves rather complicated differen- 
tial geometryS and most DNA-related theoretical studies of such models assumed circular 
cross-sections and focused on fluctuations around the straight rod configuration&il. 

Recently, we reported a study of the effect of thermal fluctuations on the statistical 
properties of filaments with arbitrary spontaneous curvature and twist@. In this work we 
present a detailed exposition of the theory and of its application to helical filaments. In 
Section |2] we introduce the description of the spatial configuration of the filament in terms 
of a triad of unit vectors oriented along the principal axes of the filament, and show that 
all the information about this configuration can be obtained from the knowledge of a set of 
generalized torsions. The elastic energy cost associated with any instantaneous configuration 
of the filament, is expressed in terms of the deviations of the generalized torsions that 
describe this configuration, from their spontaneous values in some given stress-free reference 
state. We use this energy to construct the statistical weights of the different configurations 
and show that the deviations of the generalized torsions behave as Gaussian random noises, 
whose amplitudes are inversely proportional to the bare persistence lengths that characterize 
the rigidity associated with the different deformation modes. We then derive the differential 



equations for the orientational correlation functions that can be expressed as averages of a 
rotation matrix which generates the rotation of the triad vectors as one moves along the 
contour of the filament. An expression for the persistence length in terms of one of the 
correlators is derived. In Section [3] we apply the general formalism to helical filaments and 
derive exact expressions for the correlators (see Appendix A) and for the effective persistence 
length of an untwisted helix. We show that the persistence length is, in general a non- 
monotonic function of the amplitudes of thermal fluctuations. We also show that in the weak 
fluctuation regime, our exact expressions for the correlators can be derived from a simplified 
long-wavelength description of the helix, which is equivalent to the incompressible rod-like 
chain mo defi and that the fluctuation spectrum is dominated by the Goldstone modes 
of this rod-like chain. Analytical expressions for the persistence length of a spontaneously 
twisted helix are derived (see Appendix B) and it is found that this length exhibits resonant- 
like dependence on the rate of twist. Finally, in Section f|we discuss our results and outline 
directions for future research. 



2. General Theory of Fluctuating Filaments 

A filament of small but finite and, in general, non-circular cross-section, is modeled as an 
inextensible but deformable physical curve parametrized by a contour length s (0 < s < L 
where L is the length of the filament). To each point s one attaches a triad of unit vectors 
{t(s)} whose component t 3 is the tangent vector to the curve at s, and the vectors ti(s) and 
t2(s) are directed along the two axes of symmetry of the cross-section. The vectors {t(s)}, 
together with the inextensibility condition dx/ds = t 3 , give a complete description of the 
space curve x(s), as well as of the rotation of the cross-section (i.e., twist) about this curve. 

The rotation of all the vectors tj of the triad as one moves from point s to point s' along 
the line, is generated by the rotation matrix R(s, s') 

=J2R ij (s,s')t j (s') (1) 
j 

The rotation matrix has the property 
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R(s,s') = R(s,s")R(s",s') (2) 
where s" is an arbitrary point on the contour of the filament. It satisfies the equation 

dR tJ (s,s') = _j2n ik (s)R kj (s,s') (3) 
c>s k 

where 

&>ij = '^ZijkUk (4) 
k 

Eijk is the antisymmetric tensor and {uik} will be referred to as generalized torsions, for 
lack of a better term. The above equations are supplemented by the "initial" condition 
Rij(s,s) —Sij, where Sy is the Kronecker delta function. The formal solution of Eq. (|j) is 
given by the "time-ordered" exponential 

R(s, s') = T s exp (- J' tfe"n(s")) = A hm + e " n(s " )As ■ ■ • e -n(«.)A. e -n( n )A, (5) 

The time-ordering operator with respect to s, T s is defined by the second equality in the 
above equation, where we broke the interval s — s' into n parts of length As each, so that 
Si = s' and s n = s'. The origin of the difficulty in calculating the above expression is that 
the matrices O(s) and Q(s') do not commute for s ^ s' (this is related to the non-Abelian 
character of the rotation group in 3d). 

Eq. (|^) is equivalent to a set of generalized Frenet equations from which one can calculate 
the spatial configuration of the filament, given a set of generalized torsions {uj k }, 

dti dt 2 dt 3 

— = u 2 t 3 - u 3 t 2 , — = -uit 3 + u 3 ti, — = U) x t 2 - uj 2 ti (6) 
ds ds ds 

Note that in the original Frenet description of space curves in terms of a unit tangent (which 
coincides with t 3 ), normal (n) and binormal (b), one considers mathematical lines for which 
it would be meaningless to define twist about the centerlineH. The Frenet equations contain 
only two parameters: the curvature k and torsion r: 



The two frames are related through rotation by an angle a about the common tangent 
direction (see Fig. 1), 

ti = boos a + nsina, t 2 = — bsina + ncosa (8) 

Substituting this relation into Eqs. (|5p and using Eqs. (|7|), we relate the generalized torsions 
{ujk} to the curvature k, torsion r and twist angle a, 

uj\ = /tcosa, u>2 = Ksma, = r + da / ds (9) 

The theory of elasticity of thin ro is based on the notion that there exists a stress-free 
reference configuration defined by the set of spontaneous (intrinsic) torsions {cuofc}. The set 
{u>ok} together with Eqs. (|3]) and (|J) (with uj k — > ^ofc) completely determines the equilibrium 
shape of the filament, in the absence of thermal fluctuations. Neglecting excluded- volume 
effects and other non-elastic interactions, it can be shown!! that the elastic energy associated 
with some actual configuration {uJk} of the filament is a quadratic form in the deviations 

Suik = — ^ok 

kT r L 

ds 

k 

where T is the temperature, k is the Boltzmann constant, and a» are bare persistence lengths 
that depend on the elastic constants and on the principal moments of inertia with respect 
to the symmetry axes of the cross-section, in a model-dependent way. Thus, assuming 
anisotropic elasticity (with elastic moduli Ei) and a particular form of the deformation, one 
obtainsii a\ = ExIi/kT, a 2 = EiI 2 /kT and 03 = E 2 (Ii + l-i)jkT where Jj are the principal 
moments of inertia. In general, the theory of elasticity of incompressible isotropic rods with 
shear modulus /i yields^! a x = 3/jIi/kT, a 2 = 3/j,I 2 /kT, and a 3 = C/kT where the torsional 
rigidity C is also proportional to fi and depends on the geometry of the cross-sectionil (for 
an elliptical cross-section with semi-axes bi and b 2 , C = TTfiblbl/(bf + b\) ). 

The elastic energy U e i ({5uk}) determines the statistical weight of the configuration {ujk}. 
The statistical average of any functional of the configuration B({uk}) is defined as the 
functional integral 



kT r L 

U el ({5u k }) = — dsY,a k 8u 2 k (10) 
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tn(i W\ JD{^ k }B({u; k })e- u ^/ kT 
{B ({ " fe})) = JD{Su k }e-^}/^ 

Calculating the corresponding Gaussian path integrals we obtain 

(6ui(s)) = 0, (Su^Su^s')) = ar%6(s - s') (12) 

We conclude that fluctuations of generalized torsions at two different points along the fila- 
ment contour are uncorrelated, and that the amplitude of fluctuations is inversely propor- 
tional to the corresponding bare persistence length. 

The statistical properties of fluctuating filaments are determined by the orientational 
correlation functions, which can be expressed as averages of the elements of the rotation 
matrix, 

faOOtjV)) = (Rij{s,J)) = J2(^k(s,s")R kj (s",s')) (13) 

k 

The last equality was written using Eq. (§), with s > s" > s'. Inspection of Eqs. @ 
and (f|), shows that R(s, s") depends only the torsions u k (si) with s > si > s", and that 
R(s", s') depends only on u)k(s2) with s" > S2> s' . Since fluctuations of the torsion in two 
non-overlapping intervals are uncorrelated (see Eq. (|I2"D), the average of the product of 
rotation matrices splits into the product of their averages: 

(Rij{s, s')) = £ (R ik {s, s")) (R kj {s", s')> (14) 

k 

In order to derive a differential equation for the averaged rotation matrix, we consider 
the limit As = s — s" — > 0. Keeping terms to first order in As we find 



d(R t3 (s,s')) 



ds k 



where the matrix A is defined as 



Y,A ik { S )(R kj {s,s')) (15) 



* / \ ,. &ik ~ (Rik{s,s - As)) 

A ik {s)= hm (16) 

As^0+ As 

Analogously to Eq. the formal solution of Eq. ( |T5| ) can be written as a time-ordered 
exponential, 



(R(s,s')> = T.exp (- jTcfc"A(s")) (17) 

In order to calculate the matrix A we expand the exponential in Eq. fl5|) to second order in 
As = s — s' and use the property of time-ordering operator 

ds 2 fl(s 1 )fl(s 2 ) + [ ds 2 n(s 2 )n(si / 

\s Jsi 



T s / dsi / ds 2 n(s!)n(s 2 ) = / ds x 

Js—As Js—As Js—As 



(18) 



In order to average this equation, we first calculate the average of the product 0(si)f2(s 2 ) 
using Eqs. §) and (O) 



(n( Sl )n(s 2 )) = (n( Sl )) (n(s 2 )) + dia g ( 7i )5( Sl - S2 ) (19) 

where diag(7j) is a diagonal matrix with elements 



Using Eqs. (18) and and keeping terms up to first order in As (upon integration, the 



contribution of (f2(si)) (f2(s 2 )) is of order (As) 2 ), yields 

Kk = iAk + ^2^iM^oi (21) 
I 

The elements of the averaged rotation matrix are simply the correlators of the triad 
vectors (see Eq. (0)). From the knowledge of the above correlators one can calculate other 
statistical properties of fluctuating filaments, the most familiar of which is the persistence 
length l p , that can be interpreted as an effective statistical segment length of a coarse- 
grained model, in which one replaces the filament by a random walk with the same contour 
length L and rms end-to-end separation (r 2 ): 

The end-to-end vector is defined as r = J Q L ts(s)ds and thus 

l p = lim \ f L ds f ds' (t 3 (s)t 3 (s')> (22) 

L^oo L JO JO 

The above equations describe the fluctuations of filaments of arbitrary shape and elastic 
properties, and in the following this general formalism is applied to helical filaments. 



3. Fluctuating Helices 

A. Untwisted Helix: Correlation Functions and Persistence Length 

Consider a helical filament without spontaneous twist, such that the generalized sponta- 
neous torsions {cuoa,} are independent of position s along the contour. In order to describe 
the stress-free configuration of such a filament, it is convenient to introduce the conventional 
Frenet triad which consists of the tangent, normal and binormal to the space curve spanned 
by the centerline, supplemented by a constant angle of twist «o which describes the orien- 
tation of the cross-section in the plane normal to the centerline. According to the general 
relation between the two frames, Eq. (|S|), loqx = n Q cosa , uj 02 = K sina and u 03 = r , 
where k and r are the constant curvature and torsion of the space curve in terms of which 
the total spontaneous curvature that defines rate of rotation of the helix about its long axis, 
is given by ujq = (k,q + Tq) 1//2 . The corresponding helical pitch is 2tit /uJq and the radius of 
the helical turn is 2kkq/uJq. We proceed to calculate the orientational correlation functions. 
Since A is a constant matrix, Eq. ([H]) yields (for S\ > s 2 ) 



In order to calculate the matrix e A< - Sl S2 - ) we first find the eigenvalues Aj of the matrix A, 
which are determined by the characteristic polynomial 




(23) 



A 3 - 



7A 2 + fj,X - v = 



(24) 



where we introduced the notations 



7 = 7i + 72 + 73 = aj 1 + a 2 1 + a 3 \ 



(25) 



(i = + 7172 + 7273 + 7i73, 

v = K o (7i cos 2 a + 72 sin 2 a ) + r 2 7 3 + 717273 



(27) 



(26) 



The solution of this cubic equation depends on the sign of the expression 
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A = 27 {v - v x f + 4 (// - 7 2 /3) 3 , M = ~7M - ^7 3 (28) 

For A < all the roots Aj are real. In this parameter range, fluctuations are strong 
enough to destroy the helical structure on all length scales. In the limit of very strong 
fluctuations when the bare persistence lengths are much smaller than the radii of curvature 
7 uq, we have Aj — > 73 and correlation functions become 

(t i (s l )t j (s 2 )) = e~^ si - s % (29) 



with S\ — S2 > 0. Eq. (p9|) shows that although angular correlations remain on length 
scales smaller than 1/A», they are identical to those of a persistent rod and do not carry any 
memory of the original helix. 

In the case A > 0, there is one real eigenvalue, Ai, and two complex ones, A2,3 = \R±icu, 
where 

Ai- ¥ -2— ^ + -, \R--j- (30) 
uj = V3 + ^P^j , K = 12 1 / 3 [9 (u - + V3A} 1/3 (31) 

It is shown in Appendix A that the diagonal orientational correlation functions take the 
form 

(ti (si) U (s 2 )) = (1 - a - c*) e- AlS + ( Ci e-^ s + cje*") e~ A * s (32) 

where s = si — S2 > 0. The complex coefficients c, are calculated in Appendix A. 
In the limit of small fluctuations, 7 <C ujq, we have 

Ai = ^ (1 - 2ci) Ti , A fi = ^c i7l , 2ci = l--|, ^ 2 = o; 2 (33) 

i ^0 

In this limit, it is easy to generalize our results for the diagonal correlators and write down 

expressions for all the orientational correlation functions: 

{t i (s 1 )t j (s 2 )) = 1 + [dij ^ cos(w s)e ~2^ £ ijk sm(u s)e R 

^0 \ u o J k u o 

(34) 



where s = S1 — S2 > 0. As expected, Eq. fl&f) satisfies the condition of orthonormality of triad 
vectors tj(si)tj(si) = 5^ (this geometric condition must be satisfied for the instantaneous 
triad vectors, not only on the average). Note that in the limit of weak fluctuations the local 
helical structure is preserved on contour distances s < A^ 1 and the period of rotation of the 
helix about its axis is given by its spontaneous value, 2-kuq 1 . 

Using Eqs. (^5|)-(pq) it can be shown that when A — > 0, the total curvature of the helix 
vanishes as u ~ A 1 / 2 . Since uo is positive for A > and vanishes for A < 0, in a loose 
sense it plays the role of an order parameter associated with helical order, and the point 
A = can be interpreted as the critical point at which a continuous helix to random coil 
transition takes place. However, although the dependence of u on the various parameters 
exhibits surprisingly rich behavior, the investigation of the transition region is of limited 
physical significance. The change of the helical period from 27ru; c j~ 1 to infinity takes place in 
the "overdamped" regime where this period is larger than the persistence length (a; < 7), 
and local helical structure can no longer be defined in a statistically significant sense. An 
approximate but more physically meaningful criterion for the "melting" transition is that a 
helix of period 2ttuj~ 1 melts when the persistence length becomes of order of this period. 

We now return to Eq. (E^) for the persistence length. Using the matrix equation 



Jq°° <isexp(— As) = A 1 and taking the appropriate matrix element we find: 

l P = 2^ (35) 

k% (71 cos 2 a + 72 sin a j + (r 2 + 7172) 73 

The above expression diverges in the limit of a rigid helix 7^ — > in which fluctuations have 
a negligible effect on the helix. Non-monotonic behavior is observed for "plate-like" helices, 
with large radius to pitch ratio, k /t . When no thermal fluctuations are present (7$ — > 0), 
the effective persistence length approaches zero. Weak thermal fluctuations "inflate" the 
helix by releasing stored length (by a mechanism similar to the stretching of the "slinky" toy 
spring) and increase the persistence length. Eventually, in the limit of strong fluctuations, 
the persistence length vanishes again (as 73" : ) because of the complete randomization of the 
filament. Note that the sensitivity to the (constant) angle of twist increases with radius to 
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pitch ratio. 



In the opposite limit of "rod-like" helices k — > 0, the effective persistence length ap- 
proaches 2/73 and therefore depends on a± and only and not on To and 03 which describe 
the twist of the cross-section about the centerline. This agrees with the expectation that 
since straight inextensible rods do not have stored length, their end-to-end distance and 
persistence length are determined by random bending and torsion (writhe) fluctuations only 
and are independent of twist. 

B. Weak Fluctuations: The Rod Like Chain Model 

From the discussion in the preceding section we expect that in the presence of weak thermal 
fluctuations, the filament will maintain its helical structure locally and that fluctuations will 
only affect its large scale conformation by introducing random bending and torsion of the 
helical well as random rotation of the filament about this axis. We now rederive the 

expressions for the correlators, Eq. fl34|), using a different approach that relates our work to 
that of previous investigators^ and, in the process, leads to important insights about the 
nature of the long wavelength fluctuations that dominate the spectrum of fluctuations in 
this regime. 

Note that in the absence of thermal fluctuations, 7^ = 0, the triad vectors tj attached to 
the helix can be expressed in terms of the space-fixed orthonormal triad {e} of vectors ej, 
where e3 is oriented along the long axis of the helix and ei and e2 lie in the plane normal to 
it (Fig. 2). It is convenient to introduce the Euler angles 4>o(s) = uos, 9 = arctan(Ko/r ) 
and ojq m terms of which the relation between the two frames is given by 




(36) 



where the rotation matrix 
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^ cos 0o sin 0o ^ 



V 



sin 0o cos 0o 
1 



(37) 



describes rotation by angle 0o(s) with respect to the e 3 axis. The matrix 



^ cos 6q — sin 6 ^ 



R 2 (-0o) 



V 



(38) 



/ 



1 

sin 6*o cos6> 

gives the rotation by angle —9 with respect to the e' 2 axis (e' 2 = R 3 [0o(s)] e 2 ), and R 3 (a ) 
is a rotation by angle ao about the e 3 axis (e' 3 = R 2 (— #o)e 3 ). Note that while the space- 
fixed e was taken as a conventional right-handed triad, we chose the helix-fixed t as a 
left-handed triad. Although this choice does not affect our previous results, it does affect 
the geometric relation between the two coordinate systems and, for consistency, we replaced 
the left-handed t by the right-handed one, t R = (— ti,t 2 ,t 3 ), in Eq. ( j36f) . 

In the presence of weak thermal fluctuations, the axis of the helix slowly bends and 
rotates in space, resulting in rotation of the triad {e}. Since with each point s on the helix 
we can associate its projection 



a = Tos/uj 



(39) 



on the long axis of the helix (see Fig. 2), the rotation of the triad {e} as one moves along 
this axis is given by the generalized Frenet equations, 



dei 
da 



zu 2 e 3 — zu 3 e 2l 



de 2 
da 



-zuie 3 + w 3 ei, 



de 3 
da 



wie 2 — vj 2 ei 



(40) 



The generalized torsions, Wi (s), are Gaussian random variables determined by the conditions 



{wi (a)) = 0, (wi (a) Wj (a 1 )) = a i 1 5 ij 5 (a - a') 



(41) 



where the constants a, should be determined by the requirement that the resulting expres- 
sions for the correlators (the averages of the elements of the rotation matrix) coincide with 
these in Eq. (^4|). A calculation similar to that in the previous section yields the correlators 
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(e* (cr) ej (a')) = exp (-74 |tr - a'\) , (42) 
where, analogously to Eq. (pUD , we have 

Using Eqs. (§!]), the correlators of the original triad {t} can be expressed in terms of the 
correlators of the {e} triad. Comparing the results with Eq. (|34]) , gives 

where the equality a\ = a 2 is the consequence of symmetry under rotation in the (ei,e 2 ) 
plane. 

The correlators (|41|) can be derived from an effective free energy which describes the 
long wavelength fluctuations of the helical filament, on length scales larger than the period 
of the helix L0q . 



kT 
~2 



a\ ( vo x + zu 2 ) + a 3 w 3 (45) 



This expression coincides with the elastic energy of a rod-like chain (RLC) introduced 
by Bouchiat and Mezard0. The persistence length a\ describes the elastic response to 
bending and torsion of the effective rod-like filament. The persistence length (33 controls the 
elastic response of the RLC to twist about its axis. As a consequence of the fluctuation- 
dissipation theorem, it also determines the amplitude of fluctuations A0 of the angle 0(cr) = 
^o* 7 / 7 ^ + A0(<r), where the correlator of the random angle of rotation about the axis of the 
RLC is given by 

[A0(o-) - A<j)((T')] 2 ) = CI3 1 \a - a'\ (46) 

In Eq. ([0]) we calculated the effective persistence lengths of this model (a^) in terms 
of the bare parameters of the underlying helical filament. In reference^ where the analysis 
begins with the RLC model, these corresponding persistence lengths were introduced by 
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hand. The difference between the two models becomes important if one considers the com- 
bined application of extension and twist: while such a coupling appears trivially in models 
of stretched helical filaments@, twist has no effect on the extension in the RLC model0, 
in contradiction with experimental observations^. Our analysis underscores the fact that 
the RLC model does not give a complete description of the fluctuating helix. Rather, it 
describes long wavelength fluctuations of the "phantom" axes {e^} which, by themselves, 
contain no information about the local helical structure of the filament. In order to recover 
this information and construct the correlators of the original helix (t$ (si) tj (S2)}, one has 
to go beyond the RLC model and reconstruct the local helical geometry using the relation 
between and the helix-fixed axes tj, Eqs. (|36|). 

In deriving the expressions for the correlators (tj(s)tj(O)) in terms of the correlators of 
the RLC model, we did not take into account the possibility of fluctuations of the twist 
angle of the cross-section of the helix about its centerline, «o ~~ ► a { s ) = a o + Aa(s). From 
the fact that the resulting correlators coincide with the exact expressions, Eq. (|34"D, we 
conclude that such fluctuations do not contribute to the correlators. This surprising result 
follows from the fact that in the weak fluctuation regime, the statistical properties of the 
helix are completely determined by the low energy part of the fluctuation spectrum. Such 
long-wavelength fluctuation modes (Goldstone modes) lead to the loss of helical correlations 
on length scales larger than all the natural length scales of the helix (s > 7 _1 ^> ujq 1 ). 
These Goldstone modes are associated with spontaneously broken continuous symmetries 
and correspond to bending {vj\ and zu 2 ) and twist (zu 3 ) modes of the RLC. It is important 
to emphasize that these modes correspond to different deformations of the centerline of the 
helix and not to twist of its cross-section about this centerline. Since the elastic energy, Eq. 
(PD, depends on the spontaneous angle of twist of the helix about its centerline through the 
combinations 5&\ = cos a — k cosq;o an d Su>2 = ft sin a — Kosinao, we conclude that the 
energy is not invariant under global rotation of the cross-section about the centerline and 
that such a rotation is not a continuous symmetry of the helix. Therefore, twist fluctuations 
of the helical cross-section are not Goldstone modes and do not contribute to the correlators 
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in the weak fluctuation limit. 

Another interesting observation is that there is no contribution from compressional modes 
to the long-wavelength energy, Eq. (pp. This is surprising since the RLC is a coarse-grained 
representation of the helix and the latter may be expected to behave as a compressible object, 
with accordion-like compressional mo des§. In order to check this point, we write down the 
spatial position of a point s on the helix as 

x(s) = x(a) + <Jx(s) (47) 

where x(cr) describes the curve spanned by the long axis of the helix and, therefore, defines 
the spatial position of the point a, Eq. (j39[), on the RLC contour. The deviation <5x(s) 
describes the rotation of the locally helical filament about this axis. Since the original 
filament is incompressible, it satisfies dyi/ds = t 3 . From Eq. fl36|) we obtain an expression 
for t3 which, upon substitution into the incompressibility condition and averaging over length 
scales ^> s ^> uj^ 1 (much larger than the inverse total curvature of the helix but 

much smaller than the radii of curvature of the RLC), yields 

-±i = e sW (48) 

The fact that the long-wavelength fluctuations of the helix satisfy the above incompress- 
ibility conditions, implies that compressional fluctuations do not contribute to the long- 
wavelength correlators. The origin of this observation becomes clear if we recall that the 
energy of the helix depends on the spontaneous curvature k and torsion r and, since com- 
pressional modes change the local curvature and torsion, they have a gap in the energy 
spectrum and their energy does not vanish even in the long-wavelength limit. We conclude 
that similarly to twist fluctuations of the helical cross-section, compressional modes are not 
Goldstone modes. 

The above deliberations have profound consequences for the elastic response of the fil- 
ament to long-wavelength perturbations, such as tensile forces and moments applied to its 
ends. Using the fluctuation-dissipation theorem, we conclude that as long as the deforma- 
tion of the filament remains small (on scale uJq 1 ), these forces and moments do not induce 
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twist of the cross-section of the helix about its centerline, and that the deformation can be 
completely described by the incompressible RLC model. 

C. Effect of Spontaneous Twist 

We proceed to calculate the persistence length of a helix whose cross-section is twisted by 
an angle ato(s) = dos about the centerline (do is a constant rate of twist). It is convenient 
to rewrite Eq. (p2|) as: 

L = lim \ j L ds' f L S ds (t 3 (s + s')t 3 (s')> (49) 

L^oo L JO JO 

Recall that the correlator in the integrand of Eq. fl49|) is simply the 33 element of the averaged 
rotation matrix, and is therefore the solution of equation Eq. fljjD , the coefficients of which 
are the elements of the matrix A(s + s') defined in Eq. ([H]). The diagonal elements of this 
matrix are constants while the non-diagonal elements are given by the expressions 

Ai 2 (s + s') = -A 2 i(s + s') = r + d , 

A 3 i(s + s') = — A 13 (s + s') = k q sin (d s + a ) , 

A 23 (s + s') = -A 32 (s + s') = k cos (a s + a ) , (50) 

where all the dependence on s' is contained in a = a (s'). 

The correlator in Eq. (|49"D decays exponentially fast with s, and thus the upper limit on 
the integral over s can be extended to infinity. Since the correlator is a periodic function of 
Q!q, the integration over s' can be replaced by that over ao and we obtain 

= p dao /- dg (t3(s)ts(s _ (51) 

JO 7T JO 

In deriving the above expression we assumed that the limit L — > oo is taken and that the 
total angle of twist is always large, Ldo 3> 2n (i.e., the product Ldo remains finite for 
arbitrarily small do). This assumption will be used in the following analysis. 

We first consider some limiting cases in which analytical results can be derived. In the 
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limit of vanishing twist rates, do ~~ ► 0, the persistence length is obtained by averaging Eq. 
(|35|) with respect to a . This yields: 

, = , m±^i (52, 

where 

7± = ( 7l ± 72 )/2 (53) 



with 7! and 7 2 defined in (|20|) . 

In the limit of large twist rates, do ~~ > 00 ? we can replace the denominator of Eq. (|35]) 
by its average with respect to ckq. This yields 

Ko7+ + (t 2 + 7i-7-)7s 
Finally, when 7l = 72 (ai = a 2 ), the persistence length becomes independent of twist and 
can be derived from either of Eqs. (|52| ) and (0), by substituting 7_ = 0. 

We now consider the case of arbitrary twist rates and fluctuation amplitudes. The 
calculation involves the solution of linear differential equations with periodic coefficients 
and details are given in Appendix B. We obtain: 

l, = r (55) 

An analytical expression for the complex function S(d ) is given in Appendix B. 

In Fig. 3 we present a three-dimensional plot of the persistence length given in units of 
the helical pitch /* = IuJq/2tit , as a function of the dimensionless rate of twist w = 2uq 1 cxq 
and of the logarithm of the bare persistence length a±, for a "plate-like" helix with large 
radius to pitch ratio k /t . Inspection of Fig. 3 shows that in the case of a circular cross- 
section with a\ = a 2 = 1000, the persistence length becomes independent of twist. With 
increasing asymmetry, a\ < a 2 , a maximum appears at vanishing twist rates, accompanied 
by two minima at do = icuo/2. The geometrical significance of the locations (do = 0, icuo/2) 
of these resonances is underscored by the observation that in the limit of vanishing pitch, a 
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ribbon-like untwisted (do = 0) helix degenerates into a ring. For a = ±u /2, the cross- 
section of a twisted helix rotates by ±7r with each period, and in the above limit the helix 
degenerates into a Mobius ring. As asymmetry increases (a\ <C 02), each extremum splits 
into a minimum and a maximum and eventually one obtains a dip at do = 0, accompanied 
by two symmetrical peaks at do — ±Uq/2. Note that the persistence length is a non- 
monotonic function of the amplitude of thermal fluctuations (i.e., of 1/ai): it first slowly 
increases and eventually decreases rapidly with decreasing a\. Several two-dimensional plots 
of the persistence length as a function of the rate of twist, for different combinations of the 
bare persistence lengths a, are shown in Fig. 4. The detailed behavior of the persistence 
length depends sensitively on the choice of the parameters: for example, in the limit of weak 
fluctuations three maxima are observed in Fig. 4, instead of a maximum accompanied by 
two minima in Fig. 3. In all cases, the locations of the extrema are determined by geometry 
only: d = 0, ±u /2. 

In order to demonstrate how the initial choice of the handedness of the helix breaks the 
symmetry between the effects of under and over-twist on the persistence length, in Fig. 5 we 
present a three-dimensional plot of the persistence length as a function of the dimensionless 
rate of twist w and of the inverse radius of curvature Kq, for helices with radius to pitch ratios 
of order unity and large asymmetry of the cross-section, a\ 3> a 2 . Note that for K /r < 1 
(rod-like helices), there is a single broad maximum at do = —uj /2. Then, at k /t ~ 1, a 
central peak appears at do = 0. This peak grows much faster than the do = —ujq/2 peak, 
with increasing kq/tq. At yet higher values of ko/t another peak appears at do = uj /2 and 
eventually the amplitudes of the two Mobius side-peaks become equal (and much smaller 
than the amplitude of the do=0 peak) in the limit of plate-like helices, k, /t ^> 1 (see curve 
1 in Fig. 4). 

What is the origin of the Mobius resonances observed in Figs. 3-5? Recall that the cal- 
culation of the persistence length of a twisted helix involves the solution of linear differential 
equations with periodic coefficients (Eqs. ( pB] ) in Appendix B). These equations were de- 
rived from linear differential equations with periodic coefficients and multiplicative random 
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noise, Eqs. (|3|) and Eqs. @, which are known to lead to stochastic resonancesBZl. Some 
physical intuition can be derived from the following argument. While the persistence length 
is a property of the space curve described by the Frenet triad, the microscopic Brownian 
motion of the filament arises as the result of random forces that act on its cross-section and 
therefore are given in the frame associated with the principal axes of the filament. Since the 
two frames are related by a rotation of the cross-section by an angle «o(s), the random force 
in the Frenet frame is modulated by linear combinations of sinao(s) and cosao(s). This 
gives a deterministic contribution to the persistence length which, to lowest order in the 
force, is proportional to the mean square amplitude of the random force and therefore varies 
sinusoidally with ±2a ( s )- The Mobius resonances occur whenever the total curvature of 
the helix ujq coincides with the rate of variation of this deterministic contribution of the 
random force, ±2d - 

4. Discussion 

In this work we studied the statistical mechanics of thermally fluctuating elastic filaments 
with arbitrary spontaneous curvature and twist. We constructed the equations for the 
orientational correlation functions and for the persistence length of such filaments. We 
would like to stress that our theory describes arbitrarily large deviations of a long filament 
from its equilibrium shape; the only limitation is that fluctuations are small on microscopic 
length scales, of the order of the thickness of the filament. Furthermore, since the equilibrium 
shape and the fluctuations of the filaments are completely described by the set of spontaneous 
torsions {u ok } and its fluctuations {Suk} respectively, our theory is set up in the language of 
intrinsic geometry of the space curves. All the interesting statistical information is contained 
in the correlators of the triad vectors {t} which can be expressed in terms of the known 
correlators of the fluctuations {Sulk}, using the Frenet equations. Since these equations 
describe pure rotation of the triad vectors, this has the advantage that fluctuations of the 
torsions introduce only random rotations of the vectors of the triad, and preserve their unit 
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norm. The use of intrinsic geometry automatically ensures that the inextensibility constraint 
is not violated in the process of thermal fluctuations and therefore does not even have to 
be considered explicitly in our approach. We would like to remind the readers that the 
formidable mathematical difficulties associated with attempts to introduce this constraint, 
have hindered the development of persistent chain type models in the past and led to the 
introduction of the mean spherical approximation in which the constraint is enforced only 
on the average, and to perturbative expansions about the straight rod limit. 

The general formalism was then applied to helical filaments both with and without twist 
of the cross-section about the centerline. In the latter case we found that weak thermal 
fluctuations are dominated by long wavelength Goldstone modes that correspond to bending 
and twist of the coarse-grained filament (the rod-like chain). Such fluctuations distort the 
helix on length scales much larger than its natural period but do not affect its local structure 
and, in particular, do not change the angle of twist of the cross-section about the centerline. 
Strong thermal fluctuations lead to melting of the helix, accompanied by complete loss of 
local helical structure. Depending on the parameters of the helix, the persistence length is 
a non-monotonic function of the strength of thermal fluctuations, and may first increase 
and then decrease as the amplitude of fluctuations is increased. Resonant peaks and dips 
in plots of the persistence length versus the spontaneous rate of twist are observed both 
for small twist rates and for rates equal to half the total curvature of the helix, phenomena 
which bear some formal similarity to stochastic resonances. 

There are several possible directions in which the present work can be extended. We 
did not consider here the effects of excluded volume and other non-elastic interactions, on 
the statistical properties of fluctuating filaments. Such an analysis requires the introduction 
of a field theoretical description of the filamentsii. While this approach is interesting in 
its own right, we expect that the excluded volume exponent for the scaling of the end- 
to-end distance of a single filament will be identical to that of a Gaussian polymer chain 
(self-avoiding random walk). However, new effects related to liquid crystalline ordering are 
expected in dense phases of such filaments. Another possible extension of the model relates 
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to the elasticity of random heteropolymers, with quenched distribution of elastic constants 
and/or spontaneous torsions!!. 

A natural application of our theory involves the modeling of mechanical properties and 
conformational statistics of chiral biomolecules such as DNA and RNA. The advantage of 
our theory is that it allows us to take into account, in an exact manner, the effects of thermal 
fluctuations on the persistence length and other elastic parameters of the filament. Thus, 
the generalization of the theory to include the effect of tensile forces and torques applied to 
the ends of the filament, is expected to lead to new predictions for mechanical stretching 
experiments in the intermediate deformation regime, for tensile forces that affect the global 
but not the local (on length scales < l p ) conformation of the filament. Measurements of the 
effect of elongation on thermal fluctuations of the molecule, can give information about its 
elastic constants, and help resolve long-standing questions regarding the natural curvature 
of dnaEE. It is interesting to compare our expression for the persistence length to that 
introduced by Trifonov et al.S who proposed that the apparent persistence length l a of DNA 
depends not only on the rigidity (dynamic persistence length 1^), but also on the intrinsic 
curvature of the molecule (static persistence length l s ). The apparent persistence length is 
given in terms of the two others as 

111 

t = t + t (56) 

'a <>s 

Note that the philosophy of the above approach is very similar to ours - we begin with 
filaments which have some given intrinsic length (spontaneous radius of curvature/torsion), 
and find that the interplay between this length and thermal fluctuations gives rise to a 
persistence length l p . In fact, taking for simplicity the case of a circular cross-section, 
ai = a 2 , our expression Eq. fl35|), can be recast into the form of Eq. (R3), with 



L = l P , Id = 2ai, l s = k 2 (71 + r 2 /7i) (57) 

Indeed, in our model, ai is the bare persistence length that determines the length scale on 
which the filament is deformed by thermal bending and torsion fluctuations. Our analog of 
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the static persistence length l s depends on the spontaneous bending rate k and diverges 
in the case of a straight filament (k — > 0), in which case l a —>■ Id- If we make the further 
assumption that twist rigidity is much smaller than the bending rigidity, <C ai, the 
static persistence length becomes independent of the bending rigidity and depends on both 
the spontaneous curvature and the twist rigidity. Note, however, that the resulting Kq 2 
dependence of l s differs from the originally proposed one (kq 

Another possible application of our theory involves a new way of looking into the protein 
folding problem. Usually, one assumes that the folded conformation of proteins is deter- 
mined by the interactions between the constituent amino-acids. A different approach, more 
closely related to the present work, would be to reverse the common logic: instead of try- 
ing to understand what kind of spatial structure will result for a given primary sequence 
of amino-acids, one can begin with a known equilibrium shape (native state) and attempt 
to identify the parameters of an effective filament (distributions of spontaneous torsions 
{ujQi(s)} ) which will give rise to this three-dimensional structure^. Knowledge about the 
fluctuations and the melting of proteins can then be used to determine the distribution of 
the bare persistence lengths {aj(s)}. While the question of whether such an approach can be 
successfully implemented in order to determine the relation between primary sequence and 
ternary structure remains open, our insights about the statistical properties of fluctuating 
filaments are clearly applicable to modeling of a— helices and other elements (e.g., j3— sheets) 
of secondary structure of proteins. 

Appendix A: Calculation of Correlation Functions 

We begin with the construction of the eigenvectors of the matrix A, defined by Eq. (pID , 
in the case A > (see Eq. (PB)), when there is one real eigenvalue Ai and two complex 
ones, \r± iuj. Expanding this matrix over its eigenvectors, we get 

Ay = XxUiUj + (X R + iuj) ViV* + (X R - iu) v*vj (58) 
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where the eigenvectors u, u, v, v (and the complex conjugates of the latter two, v* and v*) 
obey the ort honor mality conditions 

3 3 3 3 3 

^UiUi = ^viv* = 1, = ^ViUi = = (59) 

j=l i=l i=l i=l i=l 

Using these conditions we can exponentiate the matrix A 

u iUj e- XlS + ViV*e~^ R+iw)s + v* Vj e- {XR - iuj)s (60) 



e- As 



'.J 



Since we are interested only in the diagonal elements of this matrix, it is convenient to 
introduce the notations 

3 

Cj = ViV*, = 1 ( 61 ) 

i=l 

In addition, substituting s = in Eq. flBTiD we get 

UiUi = 1 - Ci - c * (62) 

In order to find the complex coefficients q we write down expressions for diagonal elements 
of the matrices A and A 2 

7i = (1 - ^ - cl) Ai + ^ (X R + iuj) + cl (\ R - iu) 

(63) 

(7< - ^i) 2 - + wgj = (1 - ^ - c*) A 2 + Q (A fl + ^) 2 + c* (A B - too) 2 

Looking for the solution of these equations in the form Cj = Re Cj + i Im Cj we get expressions 

for real and imaginary parts of complex parameters q 

2R = -if + gfj (Ai + Ag) + 2A fl A! + ^ 2 - c4 

w 2 + (Ai-A B ) 2 ' ( 64 ) 

2a; Im Cj = Ai — 7i + 2 (Ah — Ai) Re Cj 

Appendix B: Persistence Length of Twisted Helix 

Since the persistence length is defined by the 33 element of the averaged rotation matrix, 
we will consider the z3 component of equation Eq. (|1|) which, using Eq. (pT3|), can be 



expressed as an equation for the corresponding correlator: 

^ = - E A ^ s + s ')9i, 9i(s, s') = Us + s')t 3 (s')) (65) 
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with initial conditions g± (0, s') = g% (0, s') = and g% (0, s 1 ) = 1. The matrix A (s + s') was 
defined in Eq. fl50|). Note that since the only s-dependent parameter of the helix is the 
angle of twist, the correlators gi(s, s') depend on s' only through the parameter ao(s') = cto 
and, in order to simplify the notation, we will omit the second argument of these functions 
in the following. 

It is convenient to introduce the complex function 

f(s) = [g l (s)+tg 2 (s)}e- i ^ s+a ^ (66) 

such that / and g 3 obey the coupled equations 

as 

^ + im = ^o\{f-n (67) 
Taking a Laplace transform of these equations, 

~ poo roo 

f{p)= f (s) e-? s ds, ~g 3 (p) = / g 3 (s) e-? s ds (68) 

JO JO 

where p is, in general, a complex parameter, we get 

(p + 7+ - zr ) / (p) + iKoh (p) = - 7 -e- 2iao f (p + 2id ) , (69) 
(p + 7s) 2s (p) + [/ (p) - /* (P)] = 1 (70) 

In deriving these equations, we used the initial conditions, / (0) = and (73(C)) = 1. Sub- 
stituting (73 from Eq. ( [TCj) into (|69D, we get a closed equation for the complex function 

/: 



K 2 

(P + 7+ - iT Q ) (p + 73) + y 



/ (p) + i«o - (P) + 7- (P + 73) e- 2iQ0 f (P + 2^o) = 



(71) 



Note that the persistence length is determined by (73 (0) which can be expressed through 
/ (0) — /* (0), Eq. (|70p. The latter functions can be calculated from Eq. (|7T|), which upon 
substituting p = —2indtQ (n integer), is recast in the standard form of difference equations, 
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a n K f (-2ina ) + 2i- K f* (-2md ) + 2 1 _b n e~ 2iao f* \-2% (n - 1) a ] = (72) 
where we defined 

a n = 1 + 2 [7+ - i (r + 2na )} (73 - 2ina ) / b n = (73 - 2ina ) / (73) 

Since the persistence length is defined as the average of #3 (0) with respect to cto, it is 
convenient to introduce dimensionless functions h n as: 

K = n r ^e 2ma °f (-2md ) (74) 

JO 2,71 

We multiply Eq. (|72|) by exp (2inao) and average it with respect to a . Defining the 
parameter e = 2 / -f^/n we rewrite Eq. (|72]) in the form 

a n h n + 2i5 n0 - h*_ n + eb n h\_ n = (75) 

in which both h n and h* m enter. In order to derive closed equations for the set of {h n } only, 
we apply complex conjugation to the above equation and change n — > —n. This yields 

a -n h -n ~ 2i5 n0 ~ K + eb n h n+1 = (76) 

Substituting the equations for h*_ n and h\_ n into ( [75] ) we find 

(a n - l/a*_ n - e 2 b n b n -i/a\^ h n + 2i(l- l/a*_ n ) 5 n0 + 
2ie5 n ib n /a*_ n + eh n+1 b n /a*_ n + eh n _ x b n l a\- n = 

Let us first consider the case n ^ 0,1. Introducing new variables y n by the equality 
h n+ i = ey n h n we find 

a n - l/a*_ n - £ 2 &„&„_i/c4_ n + £ 2 y n b n /a*_ n + y~\b n /al_ n = (78) 
We now define 

A n = (a n - l/a* n ) <A_Jb n - e 2 b n „ u B n = a*_ n /a*_„ (79) 

and get the following recurrence relation, valid for n = 2, 3, ... 
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A n + l/y n -i + e 2 B n y n = (80) 

We now take n = 2 in the above equation, and solve for yi in terms of y 2 . Repeating this 
procedure (expressing y 2 in terms of y 3 , etc.) we can write the solution as a continued 
fraction 

Vl = -1/ (A 3 - e 2 B 2 / (A 3 - e 2 B 3 / {A 4 - ■ ■ ■))) (81) 



Now consider the case n — 1 in Eq. (|77f). Using the definitions of A\ and Bi, Eq. ([79|). 
it can be recast into the form: 

(A x + e 2 B iyi ) h x + 2ie + eh = (82) 
In order to obtain a closed equation for h , we return to Eq. (|76|) with n = 0, 
a*/i* - 2i - /i + £& /ii = (83) 
Eliminating hi from the above two equations we find 

h = -2i + Eh* (84) 
where, using Eq. flSTf), 5 can be represented as a continued fraction: 



a*/(l + £V(^i+£ 2 £i?/i)) (85) 
ay (l + £ 2 6 / - e 2 Bi/ (A 2 - e 2 B 2 / (A 3 - ■ ■ •)))) (86) 



The solution of Eq. §§ is: 

h o = ~ 2i Y ~pp ( 87 ) 

Recall that /io was defined as the integral over a of the function / (0) (Eq. (|74|) ) which, 
in turn, determines the Laplace transform at p = of the correlator g 3 that appears in the 
definition of the persistence length, Eq. (|5l|) . Collecting the above expressions we find: 

l P = C ^rh (0) = ^ [l - 4 " d = — ^ rpi (88) 

JO 7T 73 L 2 J 1 + (» - 1) + (-* - 1) 
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Figure captions 

Figure 1: Schematic drawing of a twisted ribbon-like filament. The vectors of the physical 
(ti, t 2 ) and the Frenet (b, n) triad can be brought into coincidence through rotation by angle 
a, about the common tangent (t 3 ). 

Figure 2: Schematic plot of section of a ribbon-like helix. The helix-fixed coordinate 
system t at contour point s' is shown. The solid line describes the associated "rod-like 
chain" to which the coordinate system e is attached at point a on its contour. The points 
a and a' on the rod-like chain are the projections of the points s and s' respectively. 

Figure 3: Three-dimensional plot of the persistence length I* as a function of the dimen- 
sionless rate of twist w and of the bare persistence length a x (logarithmic scale), for a helical 
filament with spontaneous curvature = 1, and torsion r = 0.01 (in arbitrary units). The 
bare persistence lengths are a 2 = 1000, and a 3 = 5000. 

Figure 4: Plot of the persistence length I* as a function of the dimensionless rate of 
twist w for a helical filament with spontaneous curvature k — 1, and torsion r = 0.01 (in 
arbitrary units). The different curves correspond to different bare persistence lengths: (1) 
ai = 100, a 2 = a 3 = 5000, (2) a x = 1, a 2 = a 3 = 100, (3) ai = 0.1, a 2 = a 3 = 10, (4) 
ai = 0.01, a 2 = a 3 = 10. A magnified view of the region of small twist rates is shown in the 
insert. 

Figure 5: Three-dimensional plot of the persistence length I* as a function of the di- 
mensionless rate of twist w and of the spontaneous curvature k , for a helical filament with 
spontaneous torsion To — 1 (in arbitrary units). The bare persistence lengths are a± = 500, 
a 2 = 1 and a 3 = 500. 
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